Bound hole states in a ferromagnetic (Ga,Mn)As environment 
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A numerical technique is developed to solve the Luttinger-Kohn equation for impurity states di- 
rectly in k-space and is applied to calculate bound hole wave functions in a ferromagnetic (Ga,Mn) As 
host. The rich properties of the band structure of an arbitrarily strained, ferromagnetic zincblende 
semiconductor yields various features which have direct impact on the detailed shape of a valence 
band hole bound to an active impurity. The role of strain is discussed on the basis of explicit 
calculations of bound hole states. 

PACS numbers: 75.50.Pp, 71.30.+h, 75.30.Hx 



I. INTRODUCTION 

Over the past decade, Mn-doped III-V semiconductors 
have been investigated by many groups. So far, some 
properties of those materials are not completely under- 
stood. At high doping concentrations ferromagnetic be- 
havior is observed. The ferromagnetism can be well ex- 
plained by a mean-field Zener model 0, H, Q which ac- 
counts for the exchange of the Mn d-electrons in combina- 
tion with the usual III-V semiconductor band structure. 
On the other hand, in the very dilute limit, the material 
is non-magnetic, with single valence band holes bound 
to Mn-acccptors. Previously, this diluted doping regime 
has been modelled using k.p-theory based on the spheri- 
cal model of Baldereschi and Lipari 0, [E[ and in a tight- 
binding-approximation Q. However, the crossover be- 
tween the opposite limits of high and low doping regimes 
has hardly been investigated theoretically. There are sev- 
eral situations where this crossover regime is important. 
For example, due to disorder, spatially separated low- 
holc-dcnsity islands with localized hole states may occur 
inside a ferromagnetic host material. Moreover in pat- 
terned structures, such as e.g. the tunneling structure in 
Rcf. [7}, one must be able to describe a system where 
both, localization and ferromagnetism, are important at 
the same time. 

In this paper we investigate localized hole states in 
ferromagnetic Gai_ x Mn x As. The usual approach to de- 
scribe single impurity states at the Mn-sites starts from 
the bulk band structrue of non-magnetic GaAs and incor- 
porates the magnetic moment of an isolated Mn- impurity 
by a spin-dependent local potential @ . In the crossover 
regime, which we aim to describe here, the ferromagnetic 
order in the host material should also be accounted for. 
We incorporate this by describing the bulk Gai_ x Mn x As 
host using the Zener model of [JQ. Since the magnetic 
properties of the acceptor arc already included in the bulk 
description, we can resort to a simple Coulomb potential 
without spin-dependence to describe the Mn-impurity. 

We have developed a numerical technique to solve the 
Luttinger-Kohn equation [f| for impurity states directly 
in k-spacc and in the presence of an arbitrarily complex 



host band structure. Our technique it sufficiently fast 
to avoid the need for using the spherical approximation 
[Ij]- Therefore all features of the k.p band structure of 
Gai_ x Mn x As (e.g. band warping) are taken into account. 
In this paper, a hydrogen-like impurity in strained ferro- 
magnetic Gai_ x Mn x As is considered, but the technique 
itself is restricted neither to a Coulomb potential nor to 
a particular ferromagnetic semiconductor. 



II. THE MODEL 

The general stationary Schrodingcr equation for de- 
scribing an impurity state inside a host material has 
the form 



ff (k,M,e) + £/(x) 



E,,\r) 



(1) 



where Ef, is the binding energy of the hole bound by 
an attractive potential U (x) originating from the charge 
of the impurity. _ffo(k,M,e) describes the host material 
band structure. For our present purposes is suffices to 
use a six band k.p- model for Hq. Since it is diagonal in 
k-space, Hq can be written by the k-dependent 6x6- 
matrix 

flo(k, M, e) = H LK (k) + H so + H pd (M) + H s (e). (2) 

Below, we describe briefly each constituent of Ho; the 
interested reader can find their explicit forms in the lit- 
erature UMM- 

- The Luttinger-Kohn Hamiltonian 0] HlkO*) de- 
scribes the kinetic energy of a hole in a crystal environ- 

Ikl 2 

ment. It is proportional to — J^/lkl) ( w ith m e ff the ef- 
fective mass in a given direction) and thus yields purely 
parabolic bands. The effective mass in each direction (de- 
fined as the inverse band curvature at the T-point) is de- 
termined by the three experimentally obtained Luttinger- 
parameters 71,72,73- 

- The spin-orbit Hamiltonian Q H so is constant in 
k-space. It separates the bands with a total angular mo- 
mentum j ' = i (split-off bands) from the bands with 
3 = I by A50 = 0.34eV. Because of the relatively large 
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FIG. 1: Pure GaAs band structure of the F7 and Ts bands 
calculated with two different methods: The lines which are 
present in the whole Brillouin-zone (black-online) result from 
a tight binding calculation [lf[ while the lines which are only 
plotted near F (red-online) are obtained from a k.p-type cal- 
culation. 

splitting, the 3 = \ bands are sometimes neglected, re- 
sulting in a four-band model. We do not use this approx- 
imation and work with the full six-band model. 

- The mean-field pd-exchange Hamiltonian 0, Q 
H p d(M.) is also constant in k-space but depends on 
the magnetization orientation, ft can be written as 
H pd (M) = 6B g M ■ S where M = M/|M| is the 
mean direction of all Mn d-spins and S is the valence 
band hole-spin operator. This Hamiltonian breaks time- 
reversal symmetry and thus lifts the two-fold Kramers- 
degeneracy. Together with the spin-orbit coupling it 
completely lifts the degeneracy at T. fn this paper, we 
use B g = l5meV. 

- The strain Hamiltonian Q H s (e) depends on the 
strain tensor e. fn this paper we only investigate non- 
sheared crystals (e xy = e yz = e zx = 0). ff necessary, 
shear-strain can be treated on the same footing. Since 
high-quality Gai_ x Mn x As is always grown pseudomorph 
on a substrate with a slightly different lattice constant, 
the most important strain is growth strain. Additional 
in-plane strains arising from partial relaxation [ici| and 
symmetry breaking effects Q can be observed and are 
investigated in this paper. 

- As a basis we initially work with the total angular 
momentum eigenstates \j,nij) where rrij is the magnetic 
quantum number with respect to a quantization axis in 
the [001]-direction. We refer to this basis in the following 
as the 'jm-basis'. 

The main contribution to the wave function comes 
from the small k regime of the Brillouin-zone. For this 
reason, we use the k.p-approximation, which is exact to 
order k 2 , and thus is best suited in this domain. This 
makes it a good approach to describe holes or electrons 
in semiconductors where all important electronic proper- 
ties refer to the band edges and a high numerical accuracy 
is essential. To check the quality of the k.p-mcthod deep 



inside the Brillouin zone we compare the band structures 
of pure GaAs calculated in k.p and in tight binding [TT| . 
Figure Q] shows that the top four bands nearly coincide. 
We will show later that the contribution of the split-off 
bands to an impurity state is very small so that we do 
not expect this difference in the two methods to lead to 
a significant effect on the results. 

One central point of our argumentation is that the host 
crystal is a ferromagnct and therefore has a completely 
different band structure than a non-magnetic material. 
Since the wave function of a bound hole is strongly influ- 
enced by the underlying band structure, the knowledge of 
those differences is important to understand the results of 
our calculations. The fcrromagnetism of the host mate- 
rial is modelled by the mean-field pd-cxchangc Hamilto- 
nian TJprf(M), which couples directly to the hole spin, fn 
figure [2] we show the difference between a non-magnetic 
(e.g. GaAs) and a ferromagnetic (e.g. (Ga,Mn)As) band 
structure. While in GaAs the top subbands of the va- 
lence band (i.e. the bands with j = 3/2) are degenerate 
at T, this degeneracy is completely lifted by H p d- At T, 
H p d separates the pure spin-up and spin-down (quantized 
in M-dircction) bands by 6B g . This energy separation 
leads to extremely anisotropic effective masses for small 
k and the usual labeling of 'heavy' and 'light' hole bands 
fails due to the complicated band mixing caused by the 
mean-field pd-exchange. Therefore we simply assign a 
number to label each band. 



III. SELF-CONSISTENCY EQUATION FOR 
IMPURITY STATES 

The explicit k.p-form Q of the stationary Schrodingcr 
equation JTJ) for a valence band hole trapped by an im- 
purity potential U(r) inside a host crystal is 

Y,HS n '(k)B n ,(k)+ J d 3 k'W(k-k')B„(k') =E b B n (k) 

n' 

(3) 

where -ffo(k) is the k.p-approximated band structure 
Hamiltonian of Gai_ x Mn x As [1, HJ. U{k) is the fourier 
transform of U (r) and Ef, is the binding energy of the 
impurity-hole system. B(k) = (£?i(k), Bq(Il)) is the 
envelope wavefunction of the bound hole in k-space. 
Once B is known, the actual wavefunction in real space 
can be obtained from the Fourier transform 

1&(r) =J2 [ d3k B n (k)e lk r u n (r) (4) 

where u n (r) = Uj >m . (r) is the lattice-periodic Bloch- 
function corresponding to the \j, rrij) state. When Hq de- 
scribes only one simple s-like band with a single isotropic 
effective mass and U(r) oc = is a Coulomb potential, 
equation ([3]) reduces to the non-relativistic hydrogen 
problem transformed to k-space. We have used this prop- 
erty to test the numerical procedure by application to this 
simpler problem which allows an analytical solution. 



3 




0.1 0.2 0.3 0.4 0.5 k . 
k || [100] 



(Ga,Mn)As 







___3_^r^ — 








M || [100] 





-100 ^ 
-150 



0.1 0.2 0.3 0.4 0.5 k . a 
k || [100] 

FIG. 2: The j = 3/2 valence bands of ferro- and non-magnetic 
materials. The top figure shows the typical non-magnetic 
GaAs band structure with heavy and light holes. The bot- 
tom figure is the typical ferromagnetic band structure of 
(Ga,Mn)As. The bands are labeled by numbers. The to- 
tal spin-splitting at F is 6B g . The dotted line (red-online) 
is the energy level to which the binding energy of the bound 
hole is referred. 



We arc particularly interested in the ground-state 
wavefunction of a single valence band hole bound to a 
negatively charged Mn impurity, which we treat as a 
point charge in first approximation. Rather than using 
an unscreened Coulomb potential we start with a gen- 
eral Yukawa potential to eliminate the singularity of the 
Fourier transform of - at k = 0. In case of the Yukawa 

r 

interaction, the singularity is cut off by a screening pa- 
rameter A. Up to a coupling-constant we have 



(/(,.) = £i: « ( »,_^^ ( 5 ) 



as the impurity potential. The limit A — > yields the 
Coulomb interaction. We will investigate the Coulomb 
limit numerically by checking the convergence with re- 
spect to small but finite A. 



In order to obtain a numerical solution of equation Q 
we rewrite it in a self-consistent form 



S n (k) 



E 



(k)fl n ,(k) + / d 3 k'U(k - k')Bn(k') 



(6) 



which can be solved by iteration. After each step, the 
binding energy E b is adjusted in order to obtain a self- 
consistent pair (£*£,, B(k)) which solves the Schrodinger- 
cquation 

Note that the J d 3 \t'U(k— k')i?(k') term requires a sum 
over all k' points for each k point . This means that for N 
sampling points in k-space each self-consistency step re- 
quires 0(N 2 ) operations. We have designed a dedicated 
code which analytically integrates out parts of the angu- 
lar integral [l2| and reduces the remaining integral to a 
matrix-multiplication in order to cope with this task. 

The method described above can in principle be ap- 
plied to any system with valence band holes bound to 
impurities in a ferromagnetic host, by setting the proper 
material parameters. In the remainder of this paper, we 
investigate the specific material system Gai_ x Mn x As in 
its ferromagnetic regime (x ~ 0.05). 



Valence band holes bound to Mn-impurities have been 
investigated in several papers [!, H, 0, EH . In all of these 
papers, the Mn-impurity has been treated as an isolated 
impurity center in a non-magnetic host. Here, we inves- 
tigate the Mn-center inside a ferromagnetic host, using 
the mean-field pd-exchange Hamiltonian to describe the 
ferromagnetism of the host material. Since the fourfold 
degeneracy of the zinc-blende valence band at the T-point 
is lifted by the ferromagnetism, the band structure and 
the bound hole wave function can no longer be approxi- 
mated by the spherical model of [4(. 

This deviation of the present model from the tradi- 
tional treatment of impurity states complicates a direct 
comparison of the results of the calculations. In the 
limit of vanishing ferromagnetism (B g — > 0), however, 
the binding energy of Baldercschi and Lipari [J] is recov- 
ered, which is a check of the validity of the numerical 
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algorithm. For B g ^ we measure the binding energy 
relative to the top of the split valence band at T (see 
fig. [2|). Note however that this energy level is 3B g higher 
than the GaAs band edge due to a total splitting between 
|- and |-spin of 6B g ~ 90meV. This fact should always 
be kept in mind when comparing the present results to 
non-magnetic impurity binding energies. 

Because the B g — > limit of the binding energy equals 
the Baldereschi-Lipari (BL) binding energy [ij , the same 
criticism about the binding energy (too small as com- 
pared with experiment) as for the BL-model is applica- 
ble here. To obtain binding energies that can be directly 
compared with experiments, we can follow the same path 
as BL, i.e. including a central cell correction. We have 
not done this yet, but we expect the central cell correc- 
tion to mainly manifest itself in an overall scale factor of 
the wave function extent. The main results of this pa- 
per (i.e. the shape of the wave function) should remain 
unchanged by the central cell correction. 



IV. BOUND HOLE WAVE FUNCTIONS 

The numerical solution of Eq. ((6]) results in an en- 
velope wave function of a bound hole, given as a 6- 
componcnt vector B(k) in k-space. The basis in which 
B is represented can be chosen by convenience. During 
the self-consistency cycle, the usual jm-basis is used. For 
analyzing the result, however, a different basis is more 
convenient. In figure [3] we show the k-space impurity 
wave function inside un-straincd (Ga,Mn)As in the usual 
jm-represcntation. In part (a) the magnetization is ori- 
ented along [001] and in part (b) along [100]. Since these 
two magnetization directions are equivalent, it should be 
possible to map both results (a) and (b) onto each other. 
The use of jm-representation, however, implies the choice 
of a ^-quantization axis. This choice prevents a trivial 
mapping and makes the jm-basis unattractive for visual- 
ization purposes. 

In order to obtain a comparison that is independent of 
the quantization axis, we represent B in the k-dependent 
basis in which the band structure Hamiltonian -ffo(k) is 
diagonal [l4[. Thus, a k-dependent basis transform T(k) 
is performed, so that 



T(k)#o(k)Tt(k) 



Ai(k) 







e 6 (k), 



(7) 



where £i(k) is the dispersion of the ith. band. The band 
order is chosen such that £i(0) > £j(0) for i < j. The 
transformed wave function B(k) = T(k)B(k) is from 
now on referred to as the eigenbasis-representation of 
the envelope wave function. Visualized in the eigenbasis- 
representation (see fig. HJ), the results can easily be seen 
to obey the proper symmetry with respect to magnetiza- 
tion orientation. 
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FIG. 3: Absolute values of envelope wavefunctions in k-space 
in jm-representation for k |j [100], [010], [001] with (a) magne- 
tization along [001] and (b) magnetization along [100]. The 
results are obtained from a calculation with 50 k-space direc- 
tions and 200 points in each direction. The abscissa (fc-axis) is 
given in units of Oq 1 with ao = 0.565 nm the lattice constant 
of GaAs. 
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FIG. 4: Absolute values of envelope wavefunctions in k-space 
in eigenbasis representation. As in figure[3]the magnetization 
orientation is [001] in (a) and [100] in (b). Here the different 
lines represent different bands. These bands are enumerated 
as in figure [2] 
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FIG. 5: Fit of a hydrogen ground state wave function model 
(blue online) to the numerical data of |B(k)| 2 (red online) for 
k in magnetization direction. 



Another important reason for using the eigenbasis- 
representation of the wave function is that it makes the 
connection between the wave function and the host band 
structure explicit: The ith component of B(k) refers di- 
rectly to the ith band of the host band structure at k. 
This means that by representing B in eigenbasis one can 
directly see from which bands the wave function derives. 

From figure 3] it is evident that the main contribution 
of the bound hole wave function derives from the first 
band (see fig. [2]), which is energetically highest at the 
T-point. This is a very general feature of the present 
model and results from the large spin-splitting caused by 
the pd-exchange Hamiltonian H p d- This property of the 
impurity state is a direct implication of the treatment of 
the host material as a ferromagnet. It also implies that 
the spin of the hole is also mainly determined by band 
1, i.e. it is antiparallel to the spin of the Mn d-electrons. 
Because of band mixing, this statment cannot be exact; 
from our calculations, we estimate the bound hole spin 
s z = —0.472 for magnetization along z-direction [l5| . 

The 3D fourier transform of a spherically symmetric 
function is mathematically trivial. The transformation 
of our k-space wave function B to real space, however, 
is more challenging, due to the broken spherical sym- 
metry. In order to obtain the shape of the bound hole 
wavefunction in real space we make use of the similar- 
ity to the hydrogen problem: the ground-state Coulomb 
wave function can be transformed to k-space easily. It is 
spherically symmetric and its radial part is 



C 



(i + 4fc 2 ) 2 



(8) 



where as is the Bohr radius and C a normalization con- 
stant. To extract the extent of the wave function in a 
given direction in real space, B^ h '(k) can be fitted di- 
rectly to the numerical data of B(k) for k in this direc- 

tion. To be precise, we fit ( J2 n =i l-^nC 4 )! ) to the hy- 
drogen wavefunction. Figure O shows the high quality of 
such a fit, which is due to the use of a Coulomb potential 
and the similarity to the hydrogen problem. The Bohr 




FIG. 6: Shape of bound hole wavefunctions. The top part is 
the result of a full k.p band structure model and the bottom 
part results from a simple effective mass treatment. On the 
left side the magnetization is oriented along [001], on the right 
side the magnetization is oriented along [100] . The colors are 
chosen so that the smallest extent is blue and the largest 
extent is red. The axes are given in units of the GaAs lattice 
constant ao- 



radii can be used as an interpolation dataset to obtain 
the Bohr radius as(i) in an arbitrary real-space direc- 
tion f . We visualize the wavefunction shape by plotting 
the 3D surface |r| = ag(f). 

Figure[B]shows the resulting shape of the wave function 
in real space corresponding to figures [3] and 2] Obviously, 
the extent of the wave functions is largest perpendicular 
to the magnetization direction. The physical explana- 
tion for this can be found by an analysis of the host 
band structure, especially the direction-dependent effec- 
tive masses at T: The Bohr radius of a hydrogen wave 
function of a particle with mass m e ff bound in a Coulomb 
potential which is screened by a dielectric constant k is 
given by 



a B 



m cS e 2 



(9) 



Setting m e ff to the effective mass of the top band in a 
particular direction and n ~ 12.4 to the dielectric con- 
stant of GaAs gives a rough approximation for the shape 
of the wave function. The result of this approximation is 
shown in the lower row of figure [5] Obviously the oblate 
shape of the wave function following the magnetization 
orientation is a direct consequence of the ferro-magnctic 
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FIG. 7: (Ga,Mn)As band structure for k in [100]-direction 
and M in [001]-direction. The line labeled by m G ff-modcl 
(blue-online) refers to the parabolic top band with curvature 
determined by the effective mass at the T-point and the line 
labeled by k.p-model (red-online) refers to the top band in full 
k.p-treatment. Note that the zero point of energy is chosen 
as the reference energy level of the binding energy (see also 
fig. [2] and text near by). 



band structure of the host material. In this simple model, 
only the T-point of the band structure is contributing to 
the hole wave function. The full model shows how details 
of the band structure away from T modify the result. In 
particular, the extent perpendicular to the magnetiza- 
tion shrinks when taking into account the band struc- 
ture away from T. This can easily be understood by 
considering the difference between a parabolic band with 
curvature determined by the effective mass at T and the 
same band calculated in full k.p, as shown by figure [7] 
The approximation underestimates the effective mass in 
[100]-direction. Since as oc m~^, this underestimation 
leads to an overcstimation of the wave function extent in 
this direction, as seen in figure [5] 

The wave functions deviate strongly from the spheri- 
cally symmetric result obtained using k.p-modelling by 
Bhattacharjee et al. [1, [|| for zero magnetization. We 
also obtain a much more pronounced anisotropy in the 
wave function than observed in the tight-binding model 
of [l6[. These differences are direct consequences of the 
ferromagnetism of the host material. The fact that the 
wave-function shape follows the magnetization orienta- 
tion is due to the pd-exchange [l6| and the effect is am- 
plified by the ferromagnetic band structure. 



easy axes, i.e. (100), for hole densities typical for trans- 
port samples. A growth strain in one of these directions 
makes this direction magnetically hard [3 . Because 
of this, we restrict the following discussion to magnetiza- 
tion orientations along one of the easy axes [100] or [010], 
or in between, and choose [001] as the growth direction. 
In figure [5] we show the wave functions for tensile strain 
[l9| . for a non-strained and a compressively strained host. 



6 a. 



3 a, 



tensile 
strain 



no 
strain 



compressive 
strain 




[010] 



FIG. 8: The effect of tensile (left) and compressive (right) 
growth strain on the impurity wave functions. The non- 
strained wave function in the middle is for comparison pur- 
poses. The magnetization is along [010]. 

Obviously, a compressive strain enhances the in-plane 
extent of the wave function in the magnetization direc- 
tion while a tensile strain produces a larger extent in the 
growth direction. Again, this effect can easily be under- 
stood by investigating the dependence of the band struc- 
ture on the strain properties of the host: A large tensile 
growth strain pushes the second and third bands to lower 
energies. This leads to a larger area of the Brillouin zone 
where the first band is not disturbed by the others. In 
this limit, the first band has light-hole character (large 
curvature at T) in the growth direction and heavy-hole 
character for in-plane directions. Since as oc m~g , this 
consideration qualitatively leads to the left wave func- 
tion in figure [8j The opposite limit of large compressive 
growth strain is slightly more complex. Here the second 
band becomes energetically close to the first band and 
strongly affects it. However, the effect of the second band 
is not the same for all k-directions. In growth direction, 
it pushes the first band up so that it gets heavier and the 
wave function shrinks in growth direction. In the plane, 
instead, the first band gets nearer to the second and so 
the curvature increases, leading to an enhancement of the 
wave function extent in [100] direction. 



V. STRAIN AND SYMMETRY-BREAKING 
EFFECTS 

A. Growth strain 

(Ga,Mn)As grown epitaxially on GaAs is compres- 
sively strained in the growth-plane, leading to a reduc- 
tion of the cubic symmetry and thus a change of mag- 
netic easy axes. Calculations 0, H[ and experiments have 
shown that unstrained (Ga,Mn)As has three magnetic 



B. In-plane strain 

The treatment of growth-strain in (Ga,Mn)As is well 
known and straightforward 0, 0, 0] ■ It breaks the sym- 
metry between the [001] growth direction and each of 
the two in-plane easy axes ([100] and [010]), while the 
symmetry between [100] and [010] is retained. Experi- 
mentally [13, HH, however, symmetry is also broken be- 
tween those two in-plane directions. No first principles 
explanation of this symmetry breaking has yet been es- 
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tablishcd. Rather, the effect has been modeled success- 
fully by a phenomenological in-plane strain term: Gould 
et al. introduced a strain term of order 0.1% and 
observed a magnetization orientation dependence of the 
partial DOS, leading to a remanent tunneling anisotropic 
magnetoresistance (TAMR) effect. Later, a very large 
amplification of this effect at low temperatures was ob- 
served [22J and explained in terms of a magnetization 
orientation driven metal insulator transition (MIT) Q. 
Even without any strain, the shape of a bound hole wave 
function depends on the magnetization orientation. 

A difference between the hole wave functions for mag- 
netization alon g tw o equivalent directions (in this case 
[100] and [010] [23j). however, can only be obtained by 
the introduction of a uniaxial symmetry breaking strain 
term along one of these directions - in analogy to what 
was needed to explain the remanent TAMR effect of ref- 
erence [IflJ • This is one reason for investigating the effects 
such kinds of strain theoretically. Another reason is for 
example the strongly anisotropic in-plane strain resulting 
from partial relaxation in patterned (Ga,Mn)As samples 

[MH1H 

A key ingredient of our investigation is therefore study- 
ing the effects of in-plane strain, which we refer to as 
ej = e xx in addition to the growth strain e g = e zz ; in the 
following we always use e g = 0.005, which is a reasonable 
value for (Ga,Mn)As layers with thicknesses in the lOnm 
range @, y[. We want to argue that for such parame- 
ters a switch of the magnetization orientation between 
the two in-plane easy axes can trigger a metal insulator 
transition (MIT) in a region of low hole density. Consider 
(Ga,Mn)As on the insulating side of the MIT. The mean 
overlap between neighbouring holes bound to impurity 
centers is too small to allow the holes to occupy spatially 
extended states. Then, assume some mechanism which 
causes the bound hole wave functions to grow in size 
such that there is sufficient mean overlap to destroy the 
localization. This implies that the material has crossed 
the MIT and now has a metallic character. We show 
that one such mechanism is a magnetization reorienta- 
tion from one in-plane easy axis to the other. 

In order to support this claim, we investigate the total 
extent of the wave function in terms of the volume V 
enclosed by the surface |r| = as(f) and the extent of 
the wave function in the growth direction in terms of 
as ([001]). The latter quantity is useful for investigating 
effects of tunneling in the growth direction ([001]) such 
as those in Q- 

Figure [5] shows the magnetization direction depen- 
dence of the wave function extent, measured in terms 
of V. For ei ~ we find the usual 4- fold symmetry as 
expected. For non- vanishing e.;, a strong uniaxial contri- 
bution appears which makes the total extent for magne- 
tization parallel to very different from the total extent 
for magnetization perpendicular to e,. For | | ~ |e 3 |, the 
ratio r(e s ) = V([100])/V([010]) is up to 2. A doubling 
in size of the bound hole states should easily be able to 
drive a material which is close to the MIT through this 
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FIG. 9: Relative extent of the bound hole wave functions in 
dependence of magnetization orientation. Each line is labeled 
by a particular e^. The thick black line is for ei = 0. AV is 
given relative to the mean volume of the e t = wave function. 

transition. 



C. Critical strain 

The shape of a bound hole wave function shows an- 
other interesting feature which only appears at - from a 
bulk point of view - exotic strain and pd-exchange pa- 
rameters. At interfaces, however, where the exchange 
coupling may be reduced, or in patterned samples, where 
strong in-plane strains are present, these parameters are 
reasonable. To be specific, this feature appears whenever 
the two top bands of the ferromagnetic host crystal are 
degenerate at k = 0. This happens, for example, if we 
apply a negative in-plane strain ~ —0.005 i.e. equal in 
magnitude to the growth strain. We then observe a dis- 
continous change of the shape and character of the wave 
function as a function of Cj. 

The discontinuity can best be seen from a comparison 
of the extent of the wave function in the growth direction 
a z B . In figure [10] (d) we show the e^-dependence of the 
quantity 

a%(M || [010]) -q^M || [100]) 
a%(M || [100]) 

where a^(M || [100]) is the extent of the wave function in 
growth direction when the magnetization is along [100]. 
At small €i we find the behaviour which we would expect 
from figure [HI i.e. the extent is larger (smaller) for M 
along [100] then for M along [010] for positive (negative) 
e;. At e, ~ —0.005, however, a sharp jump as a function 
of ti appears. We call the abscissa at which this jump 
appears the critical strain. Its exact value depends on 
other parameters like the pd-exchange coupling constant 
B g or the growth strain. 

Not only the growth direction extent of a hole wave 
function is affected by this effects. The overall shape of 
the wave function changes discontinously as a function 
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FIG. 10: Bound hole wave functions near critical strain: (a) 
and (b) show the wave functions for ei smaller than criti- 
cal and larger than critical respectively with magnetization 
orientation along [100]. The wave function in (c) is for mag- 
netization along [010] and strain near its critical value, (d) 
shows the relative difference of the extent of the bound hole 
wave function in growth direction ([001]) for magnetization in 
[100] and [010] direction. 



of ti at the critical strain, provided M || [100]. Figure ITD1 
shows wave functions for near its critical value, (a) and 
(b) show the difference for both sides of the criticality. 
The extent in magnetization direction is large for above 
criticality, while the extent in the plane perpendicular to 
the magnetization is smaller. For magnetization along 
[010], however, there is no such sharp switch - (c) shows 
the shape of the wave function in this case. This means 
that for ti larger than its critical value the bound hole 
wave function no longer has an oblate shape following 
the magnetization orientation. 

Let us discuss the underlying principle of the concept 
of critical strain: As mentioned above, the bound hole 
wave function derives mainly from the first band, i.e. the 
band which is energetically highest at T. At critical strain 
and M || [100], however, the first and second band are 
degenerate. In this situation, a small variation of lifts 
the degeneracy and the wave function 'switches' to one 
of the two bands - which one depends on the sign of the 
variation. The k-dependence of these two bands is in 
general very different. In figure [TT1 the most important 
charatcristic for the bound hole wave function, i.e. the 
effective mass of the first band, is shown. Obviously, the 
shapes of the wave functions in figure [TU] (a) and (b) are 
again approximately related to the sub- and super-critical 
strain effective masses by the relation as oc m~g. 

For magnetization along [010] there is no degeneracy 
at k = of the first two bands at the critical strain. 



o 
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FIG. 11: k-direction dependence of the effective mass of the 
first band at the F-point for magnetization along [100] on the 
two sides of the strain criticality. The axes are given in units 
of the free electron mass. 



A small change in e; is therefore not able to change the 
band ordering at T and no sharp switching occurs. 



VI. CONCLUSION AND OUTLOOK 

We have investigated valence band holes bound to ac- 
tive Mn-impurities inside (Ga,Mn)As in the ferromag- 
netic regime. Our main focus was the spatial extent of 
the bound hole in dependence of magnetization direction 
and strain properties of the host material. A numeri- 
cal technique has been developed to solve the Luttinger- 
Kohn equation directly in k-spacc for an arbitrarily com- 
plex host band structure. This technique is therefore 
not limited to (Ga,Mn)As but can be applied to many 
ferromagnetic systems which are describable within the 
k.p-approximation. The resulting wave functions differ 
strongly from the spherically symmetric results obtained 
by fi||. 

One of the main results is that the wave function de- 
rives mainly from the band of the host band structure 
which has highest energy at the T-point. The qualitative 
appearance of the bound hole wave function is an oblate 
shape which follows the magnetization orientation. The 
band ordering, which can be affected by external param- 
eters like strain, is extremely important for the shape of 
the impurity state wave function. Changing the order of 
the bands at T, e.g. by strain, leads to large differences 
in the shape and in the mean extent of the wave func- 
tions. With the help of this mechanism a magnetization 
switched metal-insulator transition has been explained 

0- 

It would be interesting in the future to study how a 
central cell correction affects the wave function of the im- 
purity states and their binding energy. Since this concept 
is normally incorporated by additional Gaussian correc- 
tions to the Coulomb-potential, the partial integrals en- 
tering our numerics are still feasible. 

Experimentally, we have seen evidence for hopping- 
like transport in various ferromagnetic semiconductor de- 
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vices. With the help of the detailed information about 
the wave function extent in arbitrary directions provided 
by the present model, we are now in a position to develop 
more detailed percolation models aimed at understanding 
the magnetoresistance properties of these devices. We 
anticipate that the present theory should allow to de- 



sign novel (Ga,Mn)As devices with tailored magnetore- 
sistance functionality. 
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